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COMPUTATIONAL MOIRE COMPENSATION 

BACKGROUND OF THE INVENTION 

1. Field of Invention 

[0001] This invention relates to halftones and hyperacuity printers. 

2. Background of the Invention 

[0002] High quality printing depends on many factors, including how much 
fidelity information can be passed to the printer from a data source or image 
generator. One purpose of printer electronics is to use image processing to convert 
high fidelity, sampled image data into a high-band width, high-resolution bit stream 
for delivery to a printer. Because many printers are binary, image processing 
operations such as halftoning and thresholding are sometimes used. The printer 
electronics can use these halftoning and thresholding algorithms to convert a 
rasterized image of grayscale pixels into a rasterized bitstream of high addressability 
binary on-off signals for a laser, for instance. These halftoning or thresholding 
algorithms can be complex mathematical algorithms, implemented with real time 
math processing, or computed in advance and stored in lookup tables for fast access. 

[0003] A resampling interpolator is typically found in the electronic 
pathway between the image generator and the output of the printer. The resampling 
interpolator is customarily used to provide resolution conversion, but can also be used 
to electronically compensate for spatial misalignment of the printer. The latter is 
called electronic registration, and requires fractional positioning of the data sample 
between rows and columns of pixels in the image as the relative position of the 
scanning light beam is continuously monitored with respect to the image data. To 
obtain a resample, the current position of the light beam is measured in units of 
addressability in both an x (fastscan) direction and a y (slow scan or process) 
direction. 

[0004] To implement electronic registration fully, even the halftoning 
process must be resampled at fractional intervals with respect the halftone dot grid. 
This is a problem because the halftone structure is typically restricted to integer grid 
positions to minimize moire. Moire is an error in intensity caused by uncompensated 
overlapping of the halftone grids with themselves or the output grid structure. In order 
to resample a halftone on a fractional grid, a special type of halftoner called an 
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"irrational halftoner" can be used, where off-grid resamples are allowed. 
Unfortunately, irrational halftoning can cause an increase in undesirable moire. 

[0005] Irrational halftoners are also used to make better halftone screen 
systems, such as for a full color screen system. These screen systems, such as those 
that can produce classic "rosette" halftone structures, require the resampling of the 
halftone grid on fractional intervals because the halftone angle or frequency may be 
irrational with respect to the scanning grid. For instance, to produce a halftone angle 
of 30 degrees, a resampling grid slope involving the square root of three is 
unavoidable. Once again, this need for irrational halftoning for classic "rosette" 
halftoning structures causes an increase in undesirable moire. Reference is made to 
U.S. Patent 5,410,414, "Halftoning In A Hyperacuity Printer", which is incorporated 
herein by reference in its entirety, for a full teaching of an irrational halftoner that is 
compatible with the moire cancellation extension provided by the invention herein. 

SUMMARY OF THE INVENTION 

[0006] Conventionally, attempts to minimize or eliminate moire include 
treating the cause of moire, such as, for example, by using a sampling period that is an 
exact multiple of a halftone cell, or by modifying screen angles and/or line spacings 
of the halftone dot functions, or treating the effects of moire such as, for example, by 
using error feedback methods to compensate for the quantization error introduced by 
thresholding. However, none of these methods quantifies the moire or determines a 
value that more or less cancels out the moire. 

[0007] The applicants have determined that it is possible to correlate the 
position of moire intensity errors in a simulated print with the actual memory 
locations used to produce the errors in the irrational halftoner. The moire in the 
simulated print can be measured, and, by iteration, the corresponding intensities in the 
halftoner memory locations can be raised or lowered to reduce and/or cancel the 
overall moire. During this process, a moire profile is produced that records the final 
compensation. 

[0008] This invention provides systems and methods that determine one or 
more parameters that define the moire between a halftone screen and an output grid 
structure. 

[0009] This invention separately provides systems and methods that 
determine moire phase angle zones over the full field of an image. 
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[0010] This invention separately provides systems and methods that create 
folded-field moire phase angle zones of the irrational halftoner memory. 

[0011] This invention separately provide systems and methods that 
determine and set forth a useful relationship between halftoner memory locations and 
moire phase angle zones in a full field image. 

[0012] This invention separately provides systems and methods for 
quantifying the moire that occurs in a halftone image. 

[0013] This invention further provides systems and methods for determining 
a value of the moire compensation amplitude value that obtains a lower observable 
moire for a given expected moire amplitude when combined with the image data to be 
halftoned. 

[0014] This invention separately provides systems and methods that 
digitally halftone image data using halftone data that incorporates moire 
compensation data. 

[0015] This invention separately provides systems and methods that convert 
continuous tone image data to halftone image data using a look-up table that stores 
moire-compensation-adjusted halftone dots. 

[0016] This invention separately provides systems and methods that convert 
continuous tone image data to halftone image data using a halftone function that 
incorporates a moire compensation portion. 

[0017] Various exemplary embodiments of the systems and methods 
according to this invention reduce, and, to the extent possible, minimize or 
controllably eliminate, moire that is caused by irrational halftone dots. In various 
exemplary embodiments of the systems and methods of this invention, the frequency 
and angle of the moire is determined. In various exemplary embodiments of the 
systems and methods of this invention, a halftoner, such as, for example, a 
hyperacuity halftoner, look-up table (LUT) is used. This look-up table has a two- 
dimensional (a, b) coordinate system that defines a halftone threshold function. In 
various exemplary embodiments of the systems and methods of this invention, a 
relationship between the frequency and angle of the moire and the (a, b) coordinates 
used to address the look-up table is determined. As the halftoning system advances a 
raster spot within the (a, b) coordinate system, even at an irrational angle, the 
corresponding amplitude of the moire is determined as a function of the (a, b) 
coordinates of the raster spot. In various exemplary embodiments of the systems and 
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methods of this invention, the values determined for the moire amplitude are used to 
cancel the known moire. 

[0018] The systems and methods of this invention generates computed 
moire amplitude within a two-dimensional halftone coordinate system to cancel the 
naturally occurring moire. In various exemplary embodiments, the moire- 
compensating halftone dots are pre-computed and stored within a halftoner, such as, 
for example, a hyperacuity halftoner. In various other exemplary embodiments, the 
moire canceling halftone dots are determined in real time. This allows higher 
accuracy to be achieved. 

[0019] The production of exact irrational halftones according to the systems 
and methods of this invention allows a standard "rosette" halftone structure in a three- 
color image to be obtained. Additionally, the irrational halftoning involved 
according to the systems and methods of this invention allows warping halftone 
images for electronic registration. 

[0020] These and other features and advantages of this invention are 
described in, or are apparent from, the following detailed description of various 
exemplary embodiments of the systems and methods according to this invention. 
BRIEF DESCRIPTION OF THE DRAWINGS 

[0021] Various exemplary embodiments of the systems and methods of this 
invention will be described in detail, with reference to the following figures, wherein: 

[0022] Fig. 1 is a schematic illustration of an exemplary embodiment of a 
high addressability pixel grid; 

[0023] Fig. 2 is one exemplary embodiment of a halftoner suitably usable 
with the systems and methods of this invention; 

[0024] Fig. 3 shows a portion of a halftoner memory simulated image 
overlaid with a halftone grid; 

[0025] Fig. 4 shows the same portion of a halftoner memory simulated 
image shown in Fig. 3 overlaid with moire phase angle zone lines; 

[0026] Fig. 5 shows a halftone memory with folded field moire lines of 2n 
displacement and an exemplary halftone dot function boundary overlaid on the 
halftone memory; 

[0027] Fig. 6 shows the halftone memory of Fig. 5 with folded field moire 
lines of 2k displacement, including a number of displacement zones between some of 
the lines of 2n displacement; 



5 Xerox Docket No. D/Al 088 

[0028] Fig. 7 is a graph of characteristics of spectrum vectors for a 
particular dot screen; and 

[0029] Figs. 8 through 1 1 illustrate a flowchart outlining one exemplary 
embodiment of a process for computationally canceling moire according to this 
invention. 

DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS 
[0030] When creating images using halftoning, halftone screens involving 
irrational screen angles and frequencies are often used. For example, a set of three 
screens at 141 dots per inch and rotated exactly 30 degrees apart are irrational. These 
screens are only approximated with traditional rational halftoning techniques, which 
degrade halftone dot angles desired in offset lithography. Unfortunately, when 
attempting to obtain exact halftone angles using irrational techniques, halftone dots 
are forced to land on fractional positions in the scanning grid, producing undesired 
moire. This moire is caused by the inability of halftone construction techniques to 
take into account quantization errors inevitable from the fractional overlap. 

[0031] In one exemplary embodiment of the systems and methods of this 
invention, the moire reduction systems and methods according to this invention are 
applied to a high addressability image system, such as, for example, that disclosed in 
U.S. Patent 6,449,346, incorporated herein by reference in its entirety. As pointed out 
in the 346 patent, high addressability refers to an imaging method where the imaging 
device can position the writing spot with precision finer than the size of the writing 
spot. For instance, a typical 600 spot per inch (spi) high addressability system may 
operate with a 40 micron writing spot (1 /600 th of an inch), an addressability of 
600/inch in the direction perpendicular to the raster lines, and an addressability of 
4800/inch in the direction of the raster lines. Systems that write high addressability 
images typically regulate a laser or similar writing device using clock modulation, 
amplitude modulation, pulse width modulation, pulse position modulation or 
equivalent procedures. 

[0032] Moreover, an image may be considered to be a high addressability 
image if one or more pixels within the image have a high addressability format. A 
high addressability formatted pixel (i.e., a high addressability pixel) is a pixel 
comprising one or more high addressability pixel events or units, where each high 
addressability pixel event corresponds to a specific spatial placement of the writing 
spot with respect to the pixel and has a value that represents a property of the writing 
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spot at that specific spatial placement. In a binary high addressability pixel, for 
example, each high addressability pixel event is a single bit indicating whether the 
writing spot is "on" or "off at the corresponding spatial placement. 

[0033] Illustratively, Fig. 1 is a diagram showing schematic representation 
of a high addressability image 10. Fig. 1 further shows the size of an exemplary pixel 
12 as well as the size of a writing spot 14. As shown in Fig. 1, each pixel 12 within 
the image 10 is a high addressability pixel having an addressability resolution in one 
dimension and comprising four high addressability pixel events or units 16 (e.g., a 4- 
times-high addressable, fast scan resolution). The circle is representative of the 
scanning spot, and represents a light intensity distribution having a Gaussian profile. 

[0034] As pointed out in the 396 patent, multi-level halftones can be 
generated in a digital front end (DFE) in any of a number of ways. For example, in a 
digital front end, multi-level halftoning can be performed using commonly available 
gray screening software, such as is generally available from Adobe®. Alternatively, 
the multi-level halftoning could be implemented in a hardware halftoner (not shown) 
coupled between a digital front end and rendering processor performing a hyperacuity 
function-based screening, such as that described in U.S. Patents Nos. 5,410,414 and 
5,937,109, each incorporated herein by reference in its entirety. Yet another 
possibility for generating multi-level halftones is the use of multi-level threshold 
arrays in a software or hardware implementation. 

[0035] In one exemplary embodiment of the systems and methods according 
to this invention, an irrational halftoner may be used to convert the resampled 
information from the image into a binary map. An irrational halftoner includes a 
halftoner memory, which stores pre-generated halftone cells, which represent 
intensity information (for example, 256 values), as well as integer screen positions for 
each cell. For example, a halftone cell may be formed of 8x8 or 32x32 memory 
locations. The intensity and screen position values are used to address the halftoner 
memory. A high-addressability multi-bit output is retrieved from the halftoner 
memory and is used to drive a grayscale modulator 400. 

[0036] Halftoning printers may include a halftoner memory and a 
modulator driver. The halftoner memory is generated based on a predetermined 
halftone dot pattern and is addressed by an intensity value and an integer portion of an 
irrational screen position. The memory outputs halftoned multi-bit data 
corresponding to the addressed location. The modulator driver converts the halftoned 
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multi-bit data into corresponding signals used to drive a laser diode. As a result, a 
mark is rendered corresponding to the multi-bit value at selected fast scan and slow 
scan positions. 

[0037] In certain printers, an address accumulator generates a plurality of 
screen positions in a selected halftone frequency at a selected halftone angle. Each 
screen position includes an integer component and a fractional component. A 
resampling interpolator determines an intensity value for each screen position. The 
halftoner memory stores a predetermined halftone dot function that comprises a 
plurality of halftone dot planes. Each halftone dot plane represents a dot intensity, for 
example, 0 through 255. The halftoner memory is addressed by the integer 
component of the screen position and the determined intensity value for the screen 
position. One such halftoner is disclosed in U.S. Patent 5,937,109, which is 
incorporated herein by reference in its entirety. 

[0038] As shown in Fig. 2, image data to be rendered is obtained by a 
halftone image data generator 100 from an image data source 200. An intensity 
determiner 1 10 determines the intensity or grayscale value of each pixel by, for 
example, by interpolation. An address accumulator 1 14 traverses the image data in a 
and b directions by generating (a,b) coordinates in increments based on a selected 
increment frequency and a selected increment angle. Each address designates a 
screen position (a,b) which has an associated intensity value from the image data that 
will ultimately be rendered as a mark. The a coordinate is parallel to the light beam 
fast scan direction and the b coordinate is parallel to the light beam slow scan 
direction when the halftone angle equals zero. 

[0039] The address information corresponds to the position of the light 
beam in the a,b plane. The a and b accumulator 1 14 tracks the trajectory of the light 
beam in a fast scan direction. The trajectory, however, moves through the memory at 
a screen angle which is, for example, 30 degrees. The a and b accumulator 
increments fractionally so that irrational screen positions are generated. The screen 
position (a,b) is sent to a halftoner 200. 

[0040] The halftoner 200 includes a halftoner memory 220 which is 
addressed by a plurality of intensity values. Each intensity value is represented by 
one halftone dot plane 224. In an 8-bit intensity example, the halftoner memory 
stores 256 intensity values, i.e., 256 halftone dot planes 224. Each halftone dot plane 
224 contains a halftone dot 226 defined by a dot shape and a dot size, which are 
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determined from a pre-generated halftone dot function. The halftone dot function is 
typically generated based on the specific physical characteristics of the printing 
system, such as the mark renderer 500, for which the halftoner 200 is used. 

[0041] Each halftone dot plane 224 includes a plurality of memory locations 
228, for example, in an 8x8 or 32x32 grid, as noted above. Each memory location 
228 represents a pixel 12 which is addressed by the screen position (a,b). Each pixel 
12 is made up of, for example, 4 high addressability units 16. Each address in the 
halftoner memory 220 includes three values, one for the a dimension, one for the b 
dimension, and one for the intensity dimension. These addresses are provided by the 
a address accumulator, the b address accumulator, and the gray scale data value. The 
binary gray level values generated by the halftoner 200 are used to control the mark 
renderer 500, such as a printer. 

[0042] The moire profiler and reducer 300, is used to repeatedly sample an 
image of constant intensity simulated by halftoner 200 to reduce or eliminate the 
moire produced by the simulation process. For this purpose, moire is defined as any 
difference in the average number of addressability units measured in certain parts of 
the image called zones, compared with the average number of addressability units 
produced in all the image zones combined. The moire profiler and reducer 300 
defines "full field" moire phase zones for the simulated image, and corresponding 
"folded field" moire phase zones for the halftoner dot function stored in the halftoner 
memory 220. Moire found in the simulated image in a particular full field zone will 
be reduced or eliminated during iteration by adjusting the halftoner dot function in the 
corresponding folded field zone in halftoner memory 220. This process is repeated 
until a desirably reduced moire image is produced for this intensity. Then this process 
is repeated for all intensities, for instance, 0 through 255. The halftoner output is used 
to drive a modulator 400 such as a digital-to-analog converter or a pulse width 
modulator or a shift register. The modulator 400 generates control signals that are 
usable to drive an element of the mark renderer 500 such as, for example, a laser 
diode. 

[0043] A two dimensional full field moire intensity distribution function, 
I M e , may be obtained as: 



[0044] I ME (a,b) = cos 2n— (a -cos(6 + </>)+ b-sm{e + <f>)) 



(1) 



where: 



9 Xerox Docket No. D/A1088 

[0045] a and b are the image halftone cell address coordinates (see the '414 
patent discussed above) as obtained from the accumulators that keep track of halftone 
dot locations, including both the integer and fractional portion of the halftone dot 
number; 

[0046] r is moire frequency; 

[0047] R is the halftone frequency 

[0048] 9 is moire angle; and 

[0049] <[> is the halftone angle. 

[0050] The full field moire intensity distribution function Ime gives access to 
the memory locations in an image and permit a determination of what moire phase 
angle zone those memory locations are in for the image and in the halftoner memory 
locations corresponding to those image memory locations. Moire intensity functions 
determine, for each moire phase angle zone, the average number of high 
addressability units per halftoner memory location that are in that moire phase angle 
zone. For example, using the full field moire intensity function I M e, a determination 
may be made that, in a reference zone, such as, for example, a first zone, which is, for 
example, zone 1 of 12 zones, 2.3 out of 4 possible high addressability units are 
addressed. 

[0051] The full field moire intensity distribution function I M e not only gives 
the position of the moire, but, looking at the zero crossings, the moire intensity 
distribution function indicates the location of the moire phase angle zones. 

[0052] There are two moire intensity distribution functions that predict this 
moire intensity distribution. The previously discussed moire intensity distribution 
function, called the full-field moire intensity distribution function, is expressed in 
terms of a full-field equation. The other moire intensity distribution function, called 
the folded-field moire intensity distribution function, uses the small changes in phase 
among the sample sets, instead of the absolute full-field coordinates that are used in 
the full-field equation. The full field moire intensity function I M e is set forth in Eq. 
(1), above. 
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[0054] 



[0053] 



The folded field intensity function I F is: 
I F ~ cos (2tt (k 3 - a-k 4 - b)) 



(2) 



where: 



[0055] a and b are the fractional portion of the image halftone cell address 
coordinates used in the full field equation above; and 
[0056] k 3 and 1q are halftone harmonics. 

[0057] The halftone harmonic k 3 is associated with halftone dot grid X, 
discussed below in Table 1 . The halftone harmonic k4 is associated with halftone dot 
grid Y, discussed below in Table 1. 

[0058] Fig. 3 shows part of an image simulated by halftoner 200. In this 
exemplary embodiment of the systems and methods of this invention, a halftone grid 
having 171.00 dots per inch and a halftone angle of 45.50 degrees is superimposed. 
The simulated image shows a full field moire frequency of 13.938 cycles per inch 
defined by full field moire lines of 2k displacement. The full field moire angle is 
shown as 220.796 degrees. The Fast scan direction j has 4800 addressability units per 
inch, and the slowscan direction k has 600 scans per inch. Also illustrated is an 
exemplary halftone dot 2222 inside of halftone cell 2000. The addressability 
coordinates a and b are also shown. Incorporated U.S. Patent No. 5,410,414 discloses 
how to calculate values for a halftone look-up table. 

[0059] Fig. 4 shows the same part of the simulated image shown in Fig. 3, 
emphasizing the n=12 moire phase angle zones into which each full field moire lines 
of 27i displacement have been divided. As noted above, n can vary significantly, and 
n=12 is only one exemplary embodiment used for illustration. 

[0060] Fig. 5 shows a halftone memory with folded field moire lines of 2n 
displacement and an exemplary halftone dot function boundary overlaid on the 
halftone memory. 

[0061] Fig. 6 shows the halftone memory of Fig. 5 with folded field moire 
lines of displacement, including a number of displacement zones between some of the 
lines of 27i displacement. Only two of the plurality of folded field moire lines of 27c 
displacement shown in Fig. 5 are shown divided into n=12 moire phase angle zones to 
keep the figure from becoming too cluttered, but it is understood that all of the folded 
field moire lines of 2n displacement are divided into n=12 moire phase angle zones. 
These zones are determined in the following manner. First, one uses equations (1) 
and (2), set forth above, to subdivide the full field and the folded into corresponding 
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zones. Note that both of the aforementioned functions (expressed in terms of 
equations (1) and (2)) are derived by taking a cosine of a 2-d functional value, 
resulting in a co-sinusoidal corrugation. To make zones of equal spacing, the cosine 
corrugation is replaced with a sawtooth corrugation repeating every 2n (using the 
modulus function) just as in a cosine case. However, the sawtooth function is linear, 
and the result, when normalized over the desired number of zones, is a series of zone 
boundaries within each 2n period with linear spacing. 

[0062] For the full field function, the zone number Z M e may be found by: 

mod [ 2/r— (a-cos(0+^)+&-sin(0+0)) 1,2* 

[0063] Z ME = _LL_! Tn >-A. n (3) 

where n is the number of zones. Taking the modulus of the first argument (the first 
argument being the expression in the numerator within parenthesis before ",27C M ) 
results in a number between 0 and 2n, and multiplying by n over 2n results in a zone 
number between 0 and the number of zones. 

[0064] For the folded field function, the zone number Z F is found by: 

where n is the number of zones. Taking the modulus (in the sense of the base of a 
number system) of the first argument (the first argument being the expression in the 
numerator within parenthesis before ",27c") results in a number (remainder) between 0 
and 2n, and multiplying by n over 2n results in a zone number between 0 and the 
number of zones. Any a,b can be mapped into a particular zone in the full field and 
the folded field. 

[0065] The halftone harmonics used in the aforenoted exemplary 
embodiment of the systems and methods according to this invention were chosen 
using a spectrum vector frequency graph as shown in Fig. 7. Generating halftones, 
using for example, a raster output scanner, involves a number of vectors, including a 
scan vector, a clock vector, a halftone vector x and a halftone vector y. The vector gj 
is the product of the scan vector fi and its harmonic k Jt The vector g 2 is the sum of gi 
and the product of the clock vector f 2 and its harmonic k 2 . The vector g3 is the sum of 
g 2 and the product of halftone vector f 3 and its harmonic k 3 . The vector g4 is the sum 
of the products of scan plus clock plus halftone vectors x and y with their respective 
harmonics. By selecting the harmonics as shown, below, i.e., ki = 2; k 2 = 0; k 3 = -5; 
and k4 = -5, one obtains a plot of moire frequency space. We are interested in those 
vector points landing that land near, but not at the origin, because those landing points 
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represent have the lowest frequencies (other than the origin, which represents dc, and 
cannot be seen by a human) and represent the maximum moire which will be 
observed in an image generated by the halftoner. 

[0066] In various exemplary embodiments of the systems and methods 
according to this invention, four separate vectors are employed in determining the full 

field moire intensity function. The first vector fj, the scan vector, is: 

/ ^ \ 



[0067] 
[0068] 

[0069] 



fi = 



0 



(5) 



The second vector f 2 , the clock vector, is: 



f 2 = 



V O J 



(6) 



[0070] The third vector f3, the halftone vector in the x direction, is: 



[0071] f 3 = R 



cos ^-<t)l 



sin 



\2 



[0072] where R is the halftone frequency. 

[0073] The last vector f 4 , the halftone vector in the y direction, is: 

^COS(7l - $\ 



[0074] U = R 



sin^-c^ 



(7) 



[0075] The resultant moire intensity g n is the sum of the products of the scan 
vector f i , the clock vector f 2 and the halftone vectors f3 and U, multiplied by their 
respective harmonics (k). That is, the resultant moire intensity g n is: 

[0076] g n = (k 1 -f 1 +k 2 -f 2 +k 3 -f 3 +k 4 .f 4 ) (8) 
[0077] where gi is the vector product of the scan vector fi and its harmonics 
ki, and is defined as 

[0078] gi-CV?). (9) 
[0079] g 2 is the vector sum of gi and the vector product of the clock vector 
f 2 and its harmonics k 2 and is defined as 

[0080] g 2 = (k 1 -f, +k 2 -f 2 ) ; and (10) 
[0081] g 3 is the vector sum of gi and g 2 and the vector product of the 
halftone vector f 3 and its harmonics k3 and is defined as 
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[0082] g 3 = (k, -f, + k 2 f 2 + k 3 -f 3 ). (11) 
[0083] In other words, the resultant moire intensity g n is the sum of the 

products of scan plus clock plus halftone vectors with their respective harmonics, i.e., 

the resultant moire vector. 



[0084] Table 1 lists characteristics of the spectrum vectors for a particular 
dot screen. 



Frequency 


Angle 


Harmonic 


Description 


Graph 
Symbol 


freq(fi) = 600 


Z(f,) = 90°. 


k, = 2 


Scan Frequency 


□ 


freq(f 2 ) = 1200 


Z(f 2 ) = 0°. 


k 2 = 0 


Clock Rate 


X 


freq(f 3 )= 171 


Z(f 3 ) = 0.777°. 


k 3 = -5 


Halftone Dot Grid X 


0 


freq(f 4 )= 171 


Z(f 4 ) = 2.347°. 


k4 = -5 


Halftone Dot Grid Y 


o 



TABLE 1 



[0085] This is illustrated in Fig. 7. In this example, the resultant moire 
intensity vector g n is determined to be 

rnn^i f- 10.522^ 

[0086] g n = 

1,-9.107 J 
[0087] where: 

[0088] the moire frequency, r, is 13.938 
[0089] the moire angle, 0, is 220.796°; 
[0090] the halftone frequency, R, is 171 ; 
[0091] the halftone angle, (|), is 45.5°; and 
[0092] the halftone phase, p, is 0. 

[0093] In the aforementioned exemplary embodiment of the systems and 
methods of this invention, the harmonic values, i.e., the k values, have been found by 
trial and error to sum to a value near the origin, shown in Fig 7. If a vector point 
lands directly on the origin or too close, the moire value is low frequency, and will 
normally be too low for the eye to perceive in the printed document. If the sum (i.e., 
the vector point sum of the harmonic vectors) lands too far away from the origin, the 
frequency will normally be too high for the eye to perceive in the printed document. 
Normally, in the range of perhaps from about 0.1 cycles per inch to about 100 cycles 
per inch will the moire frequency be such that the eye will notice the moire and see it 
as objectionable, although the exact frequency range may be somewhat larger. The 
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technique for adding vectors with various k harmonics in moire frequency space is a 
well-known technique in the printing industry, and is explained in detail in various 
books on the subject of moire or halftoning. See for example, "The Theory of the 
Moire Phenomenon" by Isaac Amidror, Kluwer Academic Publishers, Chapter 2. In 
the exemplary embodiment illustrated in Fig. 7, we have found objectionable moire of 
13.938 cycles per inch by using the k values 2,0,-5,-5, as in the above table. Other 
objectionable moire frequencies and angles may be found using other harmonic (k) 
values. In other words, in a given image, the overall moire profile may be made up of 
a number of component moire profiles at different frequencies. Any and all 
objectionable component moire profiles at different frequencies and angles may be 
reduced and/or eliminated or canceled according to the systems and methods 
according to this invention, simultaneously, in the context of modifying an overall 
moire profile with a single halftoner with additional computations using a number of 
different k values in the full field and/or folded field equations set forth above. 

[0094] One exemplary embodiment of a method according to this invention 
is outlined in the flowchart of Figs. 8-11. Control starts in step S1000 and proceeds to 
step SI 010, where the moire vector sum which lands nearest the origin without 
actually hitting the origin is found, the moire sampling functions values are obtained, 
and the moire phase angles are divided into zones for a number of intensity levels in 
the halftoner 200. For example a complete 360° cycle of moire phase angles may be 
divided into 12 zones of 30° each of moire phase angle. Other parameters which are 
used in determining this moire profile are the halftone screen frequency, the halftone 
screen angle <J>, the halftone screen phase p, the number of grayscale levels gi, the 
number of addressability units in the halftone Af, the number of facets per inch or the 
number of halftone dots per inch dp, and the modulus m, which is the number of bits 
used for a and y addressing in a halftone screen cell. 

[0095] Next, control proceeds to step SI 020 where all moire profiles for all 
zones and intensities are initialized to zero. During iteration for a specific intensity, 
the n moire profiles for that intensity are updated after each loop, and are intended to 
store the final result for the best profile to cancel the target moire. 

[0096] Control then proceeds to step SI 030, where a particular image 
intensity level i is selected (where i = 0 to 255, for example). Next, in step SI 040, an 
average moire profile is determined for this intensity i. Next, in step SI 050, the next 
halftoner memory location coordinate locations a and b are selected. Next, in step 



1 5 Xerox Docket No. D/Al 088 

SI 060, a folded-field equation zone number is determined for this addressability 
location, a, b. 

[0097] Next, in step SI 070, the zone number is used to access the moire 
profile for that zone. Reference is made in this regard to Figs. 4 and 6, which clearly 
illustrate that each moire phase angle zone, such as, for example, zone n=l, occurs in 
each of the plurality of zones of 2n moire lines of displacement found in a given 
image or image portion. Then control proceeds to step SI 080 where the halftone 
density is modified or servoed in the halftoner memory using negative feedback from 
the moire profile and the average moire profile. In one exemplary embodiment, a 
negative feedback formula comprises adjusting the halftoner's density for that zone by 
a fraction of the current moire profile and the average profile. By using the average 
profile in the adjustment process, the convergence caused by the servoing is made 
more stable than it would be without the average profile. 

[0098] Next, in step SI 090, a determination is made of the halftone dot 
high addressability value using the servoed (modified) density. Then, in step SI 100, 
the high addressability values for the desired halftone dot are stored in the halftoner 
memory at address a,b. Next, in step 1 1 10, a determination is made whether values 
for all halftoner memory location addresses for the given halftoner intensity value 
have been stored. If not, control returns to step SI 050. Otherwise, control proceeds 
to step SI 120 

[0099] In step SI 120, a simulated image is generated in computer memory, 
e.g., general computer memory, based on the modified (servoed) halftoner memory 
data. That image may, optionally, be printed and viewed to determine what moire 
intensity reduction appears to have occurred. 

[0100] Next, in step SI 130, temporary moire profile registers and number- 
of-entries counter are zeroed for each of the zones. (Temporary registers may be 
used to store the profile during the looping process while making the simulated image. 
The results of the process, i.e., the halftone dot values for all high addressability 
memory locations (a,b), are transferred to the profile registers for each particular 
intensity upon completion of this loop.(steps SI 050 through step SI 1 10 -No). Then, 
in step SI 140, the next fastscan j and pixel k is obtained, and the halftone positions a 
and b for this j and k are determined in step SI 150. Next, in step SI 160, the halftoner 
memory is indexed to fetch the high addressability unit (pixel event) count for this 
value of i, a, b. Then, in step 1 170, the full field zone number is computed using the 
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fractional portion of a,b and the full field equation. Next, in step SI 180, the high 
addressability unit (pixel event) count is added to the moire intensity profile register 
for this moire phase angle zone, and a number-of-entries counter for this zone is 
incremented in step SI 190. 

[0101] Next, in step SI 200, a determination is made whether this process 
has been completed for all addressability units in the fastscan and slowscan directions. 
If not, control returns to step SI 050 to get the next scan j and/or pixel k and repeats 
the process. If so, control proceeds to step S 1210 to update new moire adjustment 
values (which constitute an inverse moire profile) for all zones using the average 
number of addressability units in each zone for the particular intensity and transfer 
them to the profile registers for this intensity. 

[0102] Next, in step SI 220, a determination is made whether the adjusted 
moire intensity profile is below a preset threshold (i.e., is acceptable) for all zones. In 
one exemplary embodiment of the systems and methods according to the invention, 
this threshold is determined by computer analysis of a simulated halftoner image 
based on the moire adjustments made using the methods described herein. In another 
exemplary embodiment, the threshold is determined visually based on a video image 
or a print image based on the moire adjustments made using the methods described 
herein. However, any known or later developed quantitative method of evaluating 
moire may be used to determine one or more moire acceptability thresholds. If not, 
control returns to step SI 040 and the process repeats. If so, control proceeds to step 
SI 230 to determine whether the moire profile adjustment process described above has 
been conducted for all intensity levels. If not, control returns to step SI 030 to select 
another intensity i from 0 to 255. If so, control proceeds to step SI 240 where the 
process ends. 

[0103] Once the moire intensity level has been compensated/reduced, there 
are a number of ways to determine the effectiveness of the moire reduction. One 
method is to rely on the halftoner memory readout, as indicated above. Another is to 
make and review one or more output images made using the moire compensation 
process. 

[01 04] While this invention has been described in conjunction with the 
exemplary embodiments outlined above, it is evident that many alternatives, 
modifications and variations will be apparent to those skilled in the art. Accordingly, the 
exemplary embodiments of the invention, as set forth above, are intended to be 
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illustrative, not limiting. Various changes may be made without departing from the spirit 
and scope of the invention. 



